SLAB PERCOLATION FOR THE ISING MODEL 
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, Abstract. For the FK representation of the Ising model, we prove that the slab 

00 ■ percolation threshold coincides with the critical temperature in any dimension 

d ^ 3. 

1. Introduction 

PLh ' Renormalization arguments are at the core of the description of microscopic sys- 

tems. In particular, they are necessary close to the critical point where perturbative 
techniques no longer apply. The rigorous implementation of renormalization tech- 
niques is model dependent and in this paper we shall focus on the g-Potts model. In 
this case, the critical temperature is characterized by a breaking of symmetry (spon- 
taneous magnetization) or by the occurrence of percolation (in the FK representa- 
tion). In principle, this characterization should suffice to obtain further information 
about the sub-critical or super-critical phases, like the control of the susceptibility, 
the classification of the phases and so on. Unfortunately, some of these issues can 
only be settled at very high or very low temperatures and otherwise they remain 
^ \ open. In fact, in the intermediate regime of temperatures, the knowledge of the ab- 

sence/occurrence of percolation might not be enough to implement a mathematical 
argument. Nevertheless, progress can be made under stronger assumptions which 
c$ , sometimes can be proven afterward. 

a 

> 



o 
o 



X 



The art of statistical mechanics is to propose the right criteria from which concrete 
results can be deduced and which should be valid in the whole sub-critical/super- 
critical regime. As an example, the Dobrushin and Shlosman strong mixing property 
[DSJ implies that the sub-critical phase is well behaved (complete analyticity ...), as 
well as the dynamics. In some instances this property can be checked: for example 
Schonmann and Shlosman |SSj have shown that it is valid in the whole uniqueness 
regime of the two-dimensional Ising model with nearest neighbor interaction. 

In this paper we address another type of hypothesis, namely the percolation in 
slabs. This concept was introduced in the context of Bernoulli percolation by Aizen- 
man, Chayes, Chayes, Frohlich, Russo [ACCFR and turned out to be a crucial tool 
to derive many facts about the geometry of open clusters (see e.g. |Glj ). Let us be 
more specific. If p c denotes the critical intensity for the Bernoulli percolation, then 
for any p > p c , the origin has a positive probability to be connected to infinity by an 
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open path in 7* d . On the other hand, one can consider the stronger constraint that 
the origin is connected to infinity in a two-dimensional slice of Z d for d ^ 3. This 
happens with positive probability for any p larger than a critical value p c . Above p c 
several procedures have been developed to analyze many physical phenomena. Thus, 
all that is needed to generalize these results to the entire super-critical regime, is to 
prove that p c and p c coincide. This conjecture was solved in the breakthrough work 
of Grimmett and Marstrand |GM| also using important ideas introduced by Barsky, 
Grimmett, Newman jBGNj . The concept of slab percolation was generalized suc- 
cessfully to the random cluster measure by Pisztora [PI] , who proposed an extremely 
powerful renormalization scheme under the hypothesis that slab percolation occurs. 
Furthermore Pisztora conjectured that for FK percolation, slab percolation should 
also be valid in the whole super-critical phase. This coarse graining provides a very 
good description of the super-critical regime and in particular, it is a crucial tool in 
the derivation of the Wulff construction EE iBl lmvT] . 

The main object of this paper is to prove that for the random cluster measure 
associated to the Ising model, the slab percolation threshold coincides with the crit- 
ical temperature. The proof uses heavily the strategy of dynamic renormalization 
introduced by Barsky, Grimmett, Newman |BGNj in the context of Bernoulli per- 
colation. The starting point of their method was the assumption of percolation in 
a half space, from which a coarse graining could be implemented. For the random 
cluster measure, the basic characterization of the super-critical regime is a positive 
probability of percolation in any finite box with wired boundary conditions. The 
main difficulty with implementing this information is created by the dependence on 
the boundary conditions. 

Our approach relies on the surface tension from which accurate estimates on 
percolation in finite size volumes can be obtained uniformly over the boundary 
conditions. This requires precise controls on the surface tension which are known 
for the Ising model, thanks to inequalities. The most important property being 
the positivity of the surface tension derived by Lebowitz and Pfister [LP] in the 
whole phase transition regime. The consequences of the surface tension estimates 
are summarized in Corollary 13 . II from which the coarse graining will be constructed. 
Besides Subsection 13.11 the renormalization procedure developed in the rest of the 
paper is valid for general random cluster measures q ^ 1. Further heuristics as well 
as the scheme of the proof will be presented in Subsection 12.31 



2. Notation and Result 

2.1. The random cluster measure. Let E d be the set of bonds, i.e. of pairs 
of nearest neighbor vertices The set Q = {0, 1} E is the state space for the 

dependent percolation measures. Given uj G fl and a bond b = G E d , we say 
that b is open if a; (ft) = 1. Two sites of Z d are said to be connected if one can be 
reached from another via a chain of open bonds. Thus, each uo G Q splits Z d into the 
disjoint union of maximal connected components, which are called the open clusters 
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of Q. Given a finite subset B C Z d we use cb(c<j) to denote the number of different 
open finite clusters of u which have a non-empty intersection with B. 

For any A C Z d we define the random cluster measure on the bond configurations 
oj G — {0, 1} Ea , where Ea is the set of bonds b G E d intersecting A. The boundary 
conditions are specified by a frozen percolation configuration tt G £l c A = Using 
the shortcut c\(u) = C\(u V tt) for the joint configuration u> V tt G E d , we define the 
finite volume random cluster measure $ A ,7r on Q\ with the boundary conditions tt 
as: 




for q ^ 1. In this paper, we will sometimes use bond dependent intensities p(6). 
When there is no risk of confusion, we drop the upper-script p and write $>\. 

The random cluster measure associated to the box A N = {—N,...,N} d will 
be denoted by Q N . The measures 3> N are FKG partially ordered with respect to 
the lexicographical order of the boundary condition tt. Thus, the extremal ones 
correspond to the free (tt = 0) and wired (tt = 1) boundary conditions and are 
denoted as <& N and <& N respectively. The corresponding infinite volume (N — > oo) 
limits $ f and $ w always exist. In the following, the same set of bonds Ea will be 
used for random cluster measures in A with free or with wired boundary conditions. 

A correspondence between the q- Potts model and the random cluster measure was 
established by Fortuin and Kasteleyn |FK| (see also |ESt IG2j ). This representation 
of the Potts model will be referred as FK representation. Of particular interest for 
us is the Ising model at inverse temperature (3 which can be related to the previous 
model by setting q = 2 and choosing the bond intensity p = 1 — exp(— 2/3). More 
precisely, the Ising model on Z d with nearest neighbor interaction is defined in terms 
of spins {(Ji}i e zd taking values ±1. Let a\ N G {±I} Ajv be the spin configuration 
restricted to Aat. The Hamiltonian associated to a\ N with boundary conditions 
ctqa n is defined by 

H((t An | a dAN ) = a i a j- a i a r 

The Gibbs measure in An at inverse temperature (3 > is defined by 
fifiX" ( aA ^ = exp ( ~ P H (?*n I a 9A N )) , 

where the partition function Z a 9An is the normalizing factor. The boundary condi- 
tions act as boundary fields, therefore more general values of the boundary conditions 
can be used. 

2.2. The slab percolation threshold. The phase transition of the random cluster 
model is characterized by the occurrence of percolation above the critical intensity 

Pc 

Vp > p c , lim (0 <-> A^) = $ p ' w (0 <-> oo) > 0. (2.2) 

N^oo 



4 



T. BODINEAU 



If (3 C denotes the critical inverse temperature of the Ising model, the FK represen- 
tation implies that p c = 1 — exp(— 2j3 c ) for q = 2. 

For d ^ 3, one may wonder if the stronger property of percolation in a slab also 
holds up to the critical value. For any integers (L, N) we define the slabs of thickness 
L as 

S L>N = {—L, L} d - 2 x {-JV, . . . , N} 2 , S L = {-L, . . . , L} d - 2 x I 2 . 
A critical value can be associated to any slab thickness L 

p c (L) = inf < p ^ 0, liminf inf $^' f (0 *-> x) > 

As the function L — > p c (L) is non- increasing, it admits a limit p c as L tends to 
infinity. The critical value p c is the slab percolation threshold and satisfies 
p c ^ p c . The main result of this paper is to prove that both critical values coincide 
ifg = 2. 

Theorem 2.1. 

• If d ^ 3 and q = 2, then p c = p c . 

• More generally for q ^ 1, if Q q is the set defined before Corollary 1 lff.il then 

]p c ,i]ce,. 

To avoid many technicalities, the Theorem has been derived for random cluster 
measures with nearest neighbor interactions. Nevertheless, the proof can be adapted 
in a straightforward manner to the case of finite range interactions. Unbounded 
range interactions would require a more delicate treatment (see e.g. |MS| IGH] ). 

As a by-product of the proof, we get 

Corollary 2.1. If d ^ 3 and q = 2, then the surface tension of the Ising model (see 
Definition Iff. 1)) is equal to at p c . 

This follows from an argument similar to the one used in |BGN| to prove the 
continuity of phase transition for half-space Bernoulli percolation. 

2.3. Heuristics and scheme of the proof. In the super-critical regime, the per- 
colation cluster can be seen as a backbone on which small open clusters are attached 
(the leaves). When the bond density approaches p c , the backbone becomes thinner 
and the structure of the leaves becomes more chaotic as the correlation length ((p) 
diverges. Nevertheless, for any p > p c , patterns with similar features are repeated 
on a scale of the order of the correlation length, and the bond configurations which 
are distant from each other by at least ((p) behave essentially independently. Thus 
the backbone of the percolation cluster has enough space to spread in any slab of 
thickness much larger than ((p). This heuristic justifies the slab percolation conjec- 
ture. 

Some information about the backbone structure is encoded in the surface tension. 
If the surface tension is positive then the probability that one face of a cube of 
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Figure 1 . A backbone with some leaves is depicted. The dashed line shows which 
bonds have to be closed to disconnect the two faces of the box. 

side length N is not connected to the opposite face decays like exp(— riV^" 1 ). This 
means that to disconnect one face from the other a number of bonds of the order 
of iV d_1 have to be closed (see Figure |2J). Thus one deduces that the intersection of 
the backbone with any hyperplan contains a density of bonds. 

The implication of the surface tension on the percolation are derived in Subsection 
13.11 For the Ising model, the positivity of the surface tension has been established 
in the whole super-critical regime |LPj . Furthermore, in the framework of the Ising 
model, several estimates on the surface tension can be derived uniformly wrt the 
boundary conditions. This independence wrt the boundary conditions will be cru- 
cial to implement the renormalization scheme. With the exception of Subsection 
13.11 the rest of the paper does not rely on the properties of the Ising model. Thus 
for any q ^ 1, the coarse graining also applies under the assumption that p belongs 
to the range of intensities Q q for which Corollary 13.11 holds. We believe that the 
estimates on the surface tension derived for the Ising model should be valid for any 
random cluster model if p > p c . 

This information about the backbone can be used to implement a dynamic renor- 
malization procedure in the spirit of the one introduced by Barsky, Grimmett, New- 
man |BGNj . Before going into the details, let us comment on the results obtained in 
the case of Bernoulli percolation. Under the assumption of half-space percolation, 
Barsky et al. proposed a renormalization scheme in two steps. First they construct 
large blocks, so that with high probability the faces of these blocks are intercon- 
nected in such a way that further connections can be launched from any faces of 
these blocks. Using the independence between disjoint blocks, they were able in a 
second step of renormalization to create an infinite open cluster by piling up these 
blocks. The simplifying feature of the half-space percolation is to decouple the dif- 
ferent blocks, the drawback being that the threshold of half-space percolation may 
not coincide with p c . Grimmett and Marstrand overcame this problem by using 
the sprinkling technique. At the cost of a small increase of bond intensity, they 
were able to control the block connections without assuming half-space percolation. 
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When considering the random cluster measure, the dependence on the boundary 
conditions adds up and it seems difficult to generalize their proof directly from the 
knowledge that percolation occurs in any finite box with wired boundary conditions. 

The surface tension estimates imply precise controls on the probability of per- 
colation in a half box uniformly wrt the boundary conditions (see Subsection 13. 3|) . 
This will be used to adapt the strategy of Barsky et al. : the main effort is devoted 
to performing the first renormalization step, i.e. to constructing coarse grained 
blocks such that with high probability they satisfy good properties uniformly wrt 
the boundary conditions. As in |BGJNj . the blocks are such that a small region on 
the underside is connected to seeds lying in every facet of the block. Nevertheless, 
the lack of independence wrt the boundary conditions has prevented us from di- 
rectly applying the ideas of |BGNj and several detours are necessary. This task is 
performed in Section 0] Once this is done the coarse grained blocks can be piled 
up as if they were independent and the second renormalization step applies as in 
BGN . For the sake of completeness, the main ideas of the geometric construction, 
including the steering and branching rules, are recalled in Section 

Finally, we would like to comment on two related works based on very different 
strategies. 

The convergence of the critical temperatures of three-dimensional slabs has been 
derived by Aizenman in the context of Ising model |A1[ IA2j . The framework is 
different from ours and in particular the slabs are of the form Z 3 x {—L,L} and 
the boundary conditions are periodic. As the proof relies on the random current 
representation of the Ising model and on reflection positivity, it seems difficult to 
relate Aizenman's results to the slab percolation threshold defined in the FK setting 
(J2.3|) and therefore to Pisztora's coarse graining. 

For very large q, the FK measure can be analyzed by means of Pirogov Sinai 
theory (see e.g. [HKMt KLMR ). At p c the transition is first order and above 
p c , one can show that the only stable phase is the one with high density of open 
bonds. Thus for large q, the slab percolation threshold becomes a trivial matter 
and renormalization techniques do not seem necessary to describe the super-critical 
phase 1 . 

3. Crossing clusters 

3.1. Surface tension estimates. In this Subsection we study the dependence on 
the boundary conditions of the surface tension. 

Let us first recall the definition of the surface tension along the coordinate axis 
e*rf. Let S be a positive parameter (typically chosen very small) and N, L be two 
integers such that N ^> L. In the following these parameters are chosen of the 
form N = 2 n ,L = 2 e ,5 = 2~ p , where n,£,p are integers. We consider the increasing 
family of rectangles 

TZ L (N, 5) = {—N, . . . , N}^ 1 x {-5N -L,...,8N + L}, 
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and simply write TZ(N, 5) when L = 0. Finally, we introduce the set 3(N, 5) of bond 
configurations for which there is no connection from the top face to the bottom face 
of TZ(N, 5) (i.e. the two faces orthogonal to e^) 

2{N, 5) = {d hot n{N, 5) j4 d top n{N,5)} . (3.1) 

Definition 3.1. The surface tension in the direction is defined by 

T ? = & ~ W=i lo s ®nN, S) (W 5)) . 

For q = 2 and p = 1 — exp(— 2/3), t/ie Jsmg counterpart is 

1 ^ ± 
r B = lim lim —- log ^^'^ , 

where Z^ NS ^ denotes the partition function with mixed boundary conditions. 

The convergence of the thermodynamic limit has been derived in jMMRj . Notice 
also that the lateral sides have a vanishing perimeter and therefore they have no 
influence on the value of the surface tension. 

For our purposes, it will be useful to consider the equivalent formulation 

r ? = as & -w* log ®£w> ( m 6) ) • (3 - 2) 

The dependence on L can be estimated by using FKG inequality and the fact that 
3(N, 5) is a non increasing event 



As the LHS and the RHS (properly renormalized) converge to the surface tension, 
the identity (jS3) holds. 



The following result enables us to compare the surface tension for different bound- 
ary conditions. The parameter L should be taken large enough in order to screen 
the boundary conditions. 

Theorem 3.1. For q — 2, then for any p £ [0, 1], 

lo S ®KHN,5) (W *)) - log $g Hm (Z(N, 5)) I < (e L + c p 5 + c d ^j N d ^ , (3.3) 

where El vanishes as L tends to infinity. 

The specificity of the Ising model will be used only for the derivation of (|3.4|) . 
Thus for any q ^ 1, we introduce the set of bond intensities Q q for which (|3.4j) holds. 
The rest of the paper will apply for any random cluster measure with p £ Q q . 

We can now state the main result of this Subsection. 
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TZ L (N,S) \K{N,6) 



K{N, S) 




L 



25N 



Figure 2. The event $(N, 6) is supported by the shaded region which is decoupled 
from H. 



Corollary 3.1. Let q = 2, then for any p > p c , there exists 5 = 5(p) and L large 
such that for any N large enough 

( {d hot n(N, 5) <-> d^U(N, S)})>\- exp(-C p iV d - 1 ) , (3.4) 

where C p is a positive constant. 

This follows from Theorem 13.11 and the positivity of surface tension which was 
proven for the Ising model by Lebowitz and Pfister |LPj when f3 is larger than (3 C . 



Proof of Theorem \3.1\ 

The estimate ()3.3|) will be obtained by interpolating the boundary conditions 
wired and free. Define the set of boundary bonds at the top and bottom faces of 
7Z L (N, 5) as 

S = {(i,i + e d ), - e d ) G E n L {N ^ s) , i,j £ TZ L (N,5); i d = -j d = L + 5N} . 

Let $fcLr N g) be the wired FK measure for which the bonds in 5 have intensity 
s instead of p. The parameter s acts as a boundary magnetic field (on the faces 
orthogonal to e d ) and interpolates between "K^l^^ (for s — p) and ®-£p/ NS \ (for 
s — 0). The latter measure is the FK measure with free boundary conditions on the 
top and bottom faces of TZ L (N, 5) and wired otherwise. 

Remark 3.1. In Subsection \2.1\ the FK measure ^l^^ was defined on the set of 

bonds intersecting 1Z L (N, 5). When s = 0, the measure &kl(ns) ^ s e 1 ua ^ t° the free 
measure conditionally to the fact that the bonds in 5 are closed. Nevertheless this 
has no impact on the probability of events which are not supported by 5 and thus the 
probability of3(N,5) is the same under ^^r NS \ or &^l(n6)' 
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log$ 



K L (N,8) 



log$ 



f.w 



K L (N,5) 



3 



(N,S)) 



(3.5) 



be 



ds 



n L (N,&) 



Z{N, 5) u b 



n JO 



Si 



(N,8) 



K L {N,8) 



At this stage, a mixing property is required. We introduce now the FK measures 
on the box B# = {— K, . . . , K}^ 1 x {0, . . . , K} with intensity s for the bonds below 
the bottom face of B^ (i.e. the bonds in — e^), id = 0,i G B^}) and wired 
boundary conditions at the bottom face. When the boundary conditions on the 
remaining faces are free the corresponding FK measure is denoted by and $^' w 
for wired. 

Proposition 3.1. Let bo be the bond (0, — e<j). For any p G [0, 1] and s G (0,p] ; 



< ejjf(s) 



(3.6) 



where Ek{s) vanishes as K diverges. 

We postpone the proof of this Lemma and first estimate (|3.5|) . As 3(N,6) is 
supported by 1Z(N, 5), it is decoupled from the bonds in S. Thus using FKG property 
one can reduce the estimates to the set 1Z L (N, 5) \ TZ(N, 5) 

log Vb im (Z(N, 5)) - log &£ (m (2(N, 5i 

ds 

)M - M 



be~ J0 



< 2N d ~ l 



s(l-s) 
p ds 



TZ L (N,S)\TZ(N,8) y 







bo, 



'(N,S)\TZ(N,5) y 

+ c d LN d ~ 2 



s(l - s) 

The final bound has been obtained by applying again FKG inequality to reduce 
further the domain of the FK measure. The rest is an upper bound of the expecta- 
tion of the terms lying at a distance smaller than L from the lateral sides of 1Z L (N, 5). 



We remark that for any s G (0,p] 
1 



s(l-s) 



si 



s(l-s) 

1 



1 — p 



s + (1 - s)q / 

Thus the dominated convergence theorem and Proposition 13.11 imply that there 
exists a sequence El vanishing to as L diverges such that 



log$ 



K L (N,6) 1 



T d-2 



(d(N, 5)) - log $g wa) (Z(N, 5)) | < 2N d - l e L + Q LiV a 

It remains to compare the FK measures ®-£l{N 8) ano - ®n L (N8) ^ m °difying the 
boundary conditions on the lateral sides. This can be achieved at a finite cost 
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proportional to the perimeter of the sides parallel to e^. The corresponding error is 
bounded by the term <5iV d_1 in the LHS of ()3.3|) . This completes Theorem 13.11 □ 



Proof of Proposition IJ.il 

The strategy is to use the Gibbssian counterpart of ()3.6|) for which a mixing prop- 
erty is known. 

Step 1. 

We denote by dM K the set of boundary vertices of M K and define the set of vertices 
lying in the bottom face of dMx as 

d s M K = {i, i d = -1, -K ^i £ ^K if t < d} . 

The bond bo links the site to the site g = (0, . . . , 0, —1). The event that is 
connected to the boundary without using the bond bo is in the wired and free case 

C w = {0 m K \ {g}} and C f = {0 d s M K \ {g}} . 

Let h = h(s) = — | log(l — s) be the positive magnetic field associated to the bond 
intensity s. 

/%'>o) = *K(0 - dM K ) = $^(l C w) + - lew) W6o ) . 

Conditioning outside {bo}, we have for any boundary condition 7r which does not 
belong to C w (see e.g. equation (3.10) in |G2j ) 

$ !5W%+(i- s ) g - 

Thus 

/*>o) = «fi;(icw) + (i-*j;(ic-)) 



s + (1 — s)g 



s + (1 - s)g 

In the same way, with free boundary conditions (i.e. with zero magnetic field) 
/%>o) = <(leO + (l-<(W 

This leads to 

^.s,w/i \ A s,f /-, \ s + (1 — s)g / h + h 

(lew) - $ Bj , (IcO = (! _ s ) g M - ^ K 

Thus we have established the following equivalence 

Lemma 3.1. Fork = — |log(l — s), if the difference (/%' + (er ) — /%' Q (co)) vanishes 
to as K diverges then there exists a sequence Ek{s) such that 



«(c w )-<(c f ) 



and limx £k(s) = 0. 
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Step 2. 

Let us prove that for any h > 

l im rfj M - /4'°( ff o) = . (3.7) 

The above thermodynamic limits are taken in the upper half space. We denote by 
fi h ' + (cr ) and fi h,0 (<Jo) the corresponding limits. 

From |MMPj . we know that the functions {°~o) , /V^^o) as well as fi h,+ (a ) 
and / u ft- ' (o"o) are analytic in the complex domain {Re(h) > \lm(h)\}. Therefore the 
functions below enjoy the same properties 

^W=/4'>o)"/4>o) and = />> ) - /< V ) . 

When h > 1, then F(h) = 0. This follows from an estimate derived by Frohlich 
and Pfister (see (2.13) and (2.21) of |FF| ) 

V/i > 1, lim (a ) - A'~(a ) = . 

From the FKG inequality, the LHS dominates Fx{h) and one obtains that F(h) = 
for any h > 1. As F is analytic it must be constant for any h > 0. Thus f|3.7|) holds. 

Step 3. 

Finally, it remains to control the expectation of u>t, Q 

*S>*) = *K(lc-"0 + - l C w)o; 6o ) 

Conditionally to C w or its complement, the distribution of o>{, is determined thus 

*s>*) = ^bT(c w ) + (i - *s;(n) s + . 

A similar relation holds in the case of free boundary conditions, this leads to 

*j>o - <w = * ( x - a -ki 1 -^ ) («( cw ) - <( cf 

Combining Lemma f3. II and ()3.7|) we conclude Proposition 13.11 □ 



3.2. Crossing clusters. Corollary 13. II can be improved by showing that a connec- 
tion from the top to the bottom of TZ L (N, S) occurs with high probability. 

Proposition 3.2. For any p G Q q , there exists 5 and L such that for any N large 
enough 

®KHN,8) ( {d hot n L (N, 5) ~ d top TZ L (N, 5)} ) > 1 - exp(-CiV d - 1 ) , (3.8) 
where C = C(p, L, 5) is a positive constant. 
Proof. 
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The parameters S and L are fixed according to Corollary 13. II We first analyze the 
connections from d top TZ(N,5) to the bottom face of TZ L (N,5) and show that there 
exists C > such that 

^(N,s) ( {d hot n L {N, 5) <-> d top TZ(N, 5)} ) > 1 - exp(-CiV d - 1 ) . (3.9) 

Given a > 0, we introduce -4^ the set of bond configurations for which there are 
at least aN 6 -' 1 sites at the level {x d = —5N + 1} connected to d top 1Z(N, 5) by open 
paths contained within {xd is — SN + 1}. There exists a > and Ci > such that 

^ W) (^) > 1 - exp(-C 1 iV d - 1 ) . (3.10) 

To see this, we first notice that the set (A%) c is supported by the bonds lying in 
{%d ^ —5N + 1}. Given any configuration lj in (A%) c , it is enough to close at most 
aN 6 -' 1 bonds in order to cut all the connections from d top TZ(N,5) to d hot, JZ(N,5). 
Therefore, there exists a constant C2 such that 

<4 W ) (C^) C ) < exp(C 2 «iV d - 1 ) & nL(m (Z(N, 5)) . (3.11) 
Using the inequality 1)3.40 . one can choose a small enough such that (J3.1U)) holds. 

The next step is to extend the connections below the level {xd = —5N +1}. For 
any bond configuration ir in A^, we denote by {x^}i the first ctiV d_1 sites (wrt the 
lexicographic order) in {xd = —5N + 1} which are connected to d top TZ(N,5). Let 
C(arW) be the event that is connected to d hot lZ L (N,5) by a straight vertical 
path, i.e. using only the edges in 

V(x®) = - {k - l)e d , - ke d )) k ^ L . 

Uniformly over the boundary conditions w outside ?(a;^), we have 

«&£ (xW) (c(s«)) > $^ (I)) (C(X«)) > c L > 0. (3.12) 

To simplify the notation, we set A = 1Z (N, 5) n {2^ ^ — <5iV}. For any bond 
configuration 7r G ^4^ supported by {x^ ^ — 5iV + 1}, we get 

<5>l({d bot n L {N,5) <9 top ft(iV,5)}) ^ $A(U C ^ (i) )) 

> 1 - *a (n (c(x w )) 

The supports of the events C(x^) are disjoint. Using repeated conditionings and 
the lower bound (|3.12|) . we get uniformly over n in A^ 

<5>l({d hot n L {N,5) 3 top ft(iV,5)} ) ^ 1 - (1 - clY^ 1 > 1 - exp(-C L aN^ 1 ) , 
for some Ci > 0. 
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Finally, combining the previous estimate and (|3.9|) . there is C > such that 
<S>nHN,6)({d hot K L (N,5) <- d^TZ(N,5)}) 

[®l({d hot K L (N,5) <-> d to »n(N,5)}) Ua(Tr)" 

^ ^hHn,s)(^n) (1 -exp(-C L aiV^ 1 )) ^ 1 -exp(-CaiV a 



Thus (J3.9)) is satisfied. The Theorem can be completed by following the same strat- 
egy to extend the connections to the top face of 1Z L (N, S). 

□ 

3.3. Half box percolation. The surface tension estimates imply not only the oc- 
currence of half space percolation, but also uniform bounds wrt the boundary con- 
ditions for the corresponding finite size events. At this stage, it will be easier to 
consider sets in the half space {xd ^ 0} instead of the rectangles TZ L (N,5). We 
define below the corresponding notation which will be used throughout the paper. 

Definition 3.2. 

(1) For any integers (£, h), a block is defined as 

M(£,h) = {-£,..., £} d ~ l x{0,...,h}. 

The interior ofM(£,h) is 

M*{£,h) = { -£+l,...,£-l} d ~ 1 x {l,...,h-l}. 

(2) The top face of the block M(£, h) will be denoted by 

T(£,h) = {-£,..., £} d ~ l x{h}. 
The top face is split into 2 d ~ 1 sub-regions 

T 2 (£,h) = {0,...,£}^ 2 x{-£,...,0}x{M, 

T 2d -i(l,h) = {-^...^y-'xih}. 

(3) Denote by S(£, h) the union of the sides ofM(£, h). It is divided into 2(d — 
l)2 d ~ 2 sub-regions 

S^h) = {£}x{0,...,£} d ~ 2 x{0,...,M, 

S 2 (£,h) = mx{0,...,£} d - 3 x{-£,...,0}x{0,...,M, 

S 2{d _ ipd -2(£,h) = H,...,0}"x{f}x{0 ) ...>}. 



For any % ^ d and x e Z d , let b l K {x) be the {d— l)-dimensional hypercube centered 
around x and orthogonal to e*j 

b K( x ) = {v e I?, yi = Xi, Vj ^ i, \yj -Xj\ ^ K] . 
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In the following we will be interested to the connections from b l K (fS) to the faces of 
M(£, h) by open paths strictly within M(£, h), i.e. paths lying in E B *(^). Recall that 
Eb*(£,/i) denotes the set of bonds intersecting M*(£, h) (see Subsection 12. 

Proposition 3.3. For any p G Q q , there exists 5 > and L such that the following 
holds uniformly in N and K ^ N 

<(N,n N) ({b d K(0) - T(N,H N )})>l-exp(-CK d - 1 ), (3.13) 
where = d~N + L and C is a positive constant. 
Proof. 

Let B N = { - N, . . . , N}^ 1 x {0} be the bottom face of M*(N, H). According 
to Proposition 13.21 the following holds for appropriate choices of 5 > and L 

&m*(n,h> n ) ( {Bn - T(N, H' N )})>1- exp(-CiV rf - 1 ) , (3.14) 

where H' N = 25N + 2L. 
We consider also 

] = {x G B N} x = 2Kj for j G Z d } . 
Then B N is covered by {J xeB {V b d K (x) and 

{B N ^T(N,H' N )}= f| {b d K (x)^T(N,H' N )} . 

These events are decreasing, therefore ()3.14|) and FKG inequality imply 

exp(-CiV d - 1 ) > J] ®i HN>H , N) ({b d K (x)^T(N,H' N )}) . (3.15) 

Applying again the FKG property, we see that for any x in B^ 

&n*(N,H' N ) t4 T(N,H' N )}) > <* (2W (K(0) T(N , H' N ) - x}) 

> ^(2N,H' N )({b d M^T(2N,H' N )}) . 

Plugging this inequality into ()3. 15|) leads to 

-C— w > log®l {2N)H , N) ({61(0) t4 T(2iV,^)}) , 



where | = (N/2K) d 1 is the cardinal of .Bj^ . This completes the Proposition 
/ 

N/2- 

□ 



with H N = Hi 
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3.4. The seeds. We introduce now the notion of seed which will be essential to 
iterate (l3~T5j) . 

Definition 3.3. 

(1) We say that b d K {x) is a seed centered at site x G Z d if all the bonds lying in 
b d K (x) are open. 

(2) For any (£, h), let Ck{£, h) be the event that 6^(0) is connected by a path of 
open bonds strictly within M(£, h) to a seed b d K (x) included in T(£, h). 

(3) For any (£,h), C l K (£,h) denotes the event that b d K {Q) is connected strictly 
within M(£,h) to a seed b^lx) included in T(£,h) and centered at site x in 
Ti(£,h), withi ^2 d -\ 

In the following, we choose p G Q q and according to Proposition 13.31 

Proposition 3.4. For any integer K there exists an integer M = M(K) and a 
height h = h(K, N) G [Hpj — M, Hn] such that uniformly over N large enough 

<(n,h n) (Ck(N, hj) > 1 - 4 exp ( - CK d ^) , 
for some positive constant C . 
Proof. 

For a given height h, define Y (£, h) as the number of sites in T(£, h) which are 
connected strictly within M(£, h) to b d K (Q) by paths of open bonds. 

Let m, M be two integers which are going to be chosen later. The sequence of 
random heights Hi is defined recursively. We set 7i = H N — M and 

H i+1 = inf {h > H h 1 ^ Y(N, h) < m} A H N . 

We first check that there exists r p = r p (m) < 1 such that uniformly in iV and M 

\fn G [2, M], $i* (N , HN) (K n < H N ) ^ r;- 1 . (3.16) 

To see this, we write for n ^ 2 

®m*(N,H N )y~Ln < H N ^j = ^2 ®n*(N,H N )(7~Ln-i = h,H n < H N ^j . 

h<H N -l 

Conditionally to the event {Ti n -i = h} any bond configuration in {H n < H N } is 
such that Y(N, h + 1) ^ 1. On the other hand, there are at most m sites of T(N, h) 
connected to b d K {fS). Thus 







^ l-c p n = r p <l 



where c™ > is a lower bond of the energetic cost for closing the bonds which could 
connect b d K {fS) to the height h + 1. By iterating the procedure, we obtain (|3.16|) . 
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For any K, there exists n = n(K, m) such that for any M large enough and 
uniformly over N, 

<*(n,h n )( E UY(N, h>m} >M-n\ ^l-2exp(-CK d ~ 1 ), (3.17) 

\h=H N -M J 

where C is the constant of Proposition 13.31 
To see this, we write 




M*(N,H N ) \ 2—< L {Y(N,h) ^m} II ] ^ ^V*(N,H N ) 
K h=H N -M 

+ ®B*(N,H N ) < H N ) . 

The first term in the LHS means that there is no connection from b d K (Q) to the 
top face of M(N,Hn) thus it is bounded by Proposition 13.31 Let n = n(K,m) be 
such that Tp~ l ^ exp(— CK d ~ l ). Then the second term is bounded by the estimate 
(Jinnj). Thus (jTrZP is satisfied. 

We need to check that there exists h G [Hn — M, Hn] such that 

^ {N ,H N )(Y(N,h)^m)^l-3ex V (-CK d - 1 ) . (3.18) 

Applying Tchebyshev inequality to the estimate ()3.17|) . we get 

H N 

(l-2exp(-CK d - 1 ))(M-n)<: ^ < H n,h n ) (Y(N,h) > m) . 

h=H N -M 

We choose M = M(K,m) such that n(K,m) < exp(—CK d ~ r )M then the above 
inequality ensures that there exists a height h in [Hn — M, Hn] such that ()3.18|) 
holds. The height h depends on K, N and m, nevertheless the parameter m will be 
fixed according to K (see (I3.19jl ). Thus, ultimately, only the dependency on K, N 
will remain. 



The estimate (J3.18J) implies that with high probability, there are at least m sites 
in T(N, h) connected to 6^(0). Thus it is enough to find at least a seed lying on top 
of one of these sites. The probability that uniformly wrt the boundary conditions 
there exists a seed of side length K is bounded from below by , with c p > 0. 
Let us choose m = m(K) such that 

(i _ c K«-y/(*Kr-> < exp( _^-i) . (3 . 19) 

Conditioning wrt the bond configurations below h which belong to the event 
{Y(N,h) m}, we can find at least m/(4i^) d_1 sites at distance larger than 4K 
from each other. Then the probability that a seed (lying in T(N, h)) is attached on 
top of one of these sites is dominated by independent Bernoulli variables. As m has 
been chosen large enough, there must be at least a seed on one of the attachment 
site with probability larger than 1 — exp(— CK d ~ x \ 

This concludes the proof of Proposition 13.41 □ 
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Proposition 13.41 has to be enhanced in order to control the position of the seeds. 

Proposition 3.5. For any large K there exists an integer M = M(K) and a height 
h = h(K, N) G [Hn — M, Hn] such that uniformly over N large enough 

Vz ^ 2 d ~\ n* { N,H N) (CM h))>l- exp ( - CK d ~i) , 
for some positive constant C. The event C Z K (N, h) was introduced in Definition ^, ft 

Proof. 

Following jRGN] . we use the symmetry of the model to write 

(<*(N,H N) ((c i K (N,h)y)y d 1 = u^ iN , HN) ((ci(N,h)r) 

< <*(N, HN )(r\( C K( N > h )) C ) = <*(N,H N )((C K (N,h)) C ) , 
3 

where we used the FKG inequality and the fact that the events C J K (N,h) are in- 
creasing. 

Using the lower bound derived in Proposition 13.41 we obtain that 

V. < 2 d -\ ®1 {n , Hn) ((CitiN, h)) c ) < 4 1 ^- 1 exp (-T^T^ 1 ) • 
This completes the proof. □ 



4. The coarse graining: first renormalization step 

Recall that p is chosen in the set Q q introduced before Corollary 13.11 

4.1. Top connections. The next step is to iterate Proposition l3.5l in order to obtain 
a connection in a box M(N, L^) for some ^ 3iV. 

Theorem 4.1. For any rj > 0, one can find K = K{rj) such that uniformly in N 
(large enough), there exists an integer M = M(K) and a height (p = <p(K,N) G 
[L N — M, L N ] such that 

where L N = L N (K) is such that 3N ^ L N - M(K) ^ <p(K, N) (see (@~U) )■ 
Proof. 

Let Hn = 8N + L, where 8 and L are chosen according to Proposition 13.31 Fix 
K large enough such that 

(1 - exp (-CA^" 1 )) 1+6/5 ^l-Ty, 

where C is the constant obtained in Proposition 13.51 Then choose N large enough 
and set h = h(K, N),M — M(K) according to Proposition 13.51 



18 



T. BODINEAU 



2N 
N 




- b d K (0) 

Figure 3. The seeds are connected according to the steering rule. 

We set 

n = 6/<5, L' N = H N + nh(K,N), L N = L' N/2 . (4.1) 

Let Ai be the event that 6^(0) is connected to a seed in T(N,ih) by an open path 
strictly contained within M(2N,ih) 

Mi < ra, Ai = {b d K (0) <-> seed cT(N,ih)} . 

We stress the fact that the event Ai differs from Ck(2N, ih) since the attachment 
site of a seed must be at distance less than N from the axis {x\ = 0, . . . , Xd-i = 0}, 
nevertheless the open path leading to this seed is allowed to use any bonds within 
M(2N,ih). K(p'(K,N) = (n+l)h(K,N) then A n+ i C C K (2N, <//). * 

Conditioning wrt rr n the bond configuration below the height nh(K, N), we get 

®B*(2N,L' N ) (Ai+l) ^ ®B*(2N,L' N ) (Ai H At+l) 

^ *B W W) (A (Wy^) " Seed C ( n + W})) > 

where ?/ n ) is defined for each 7r n as the attachment site of a seed in T(iV, nh) (if 
there are many such sites choose the smallest one wrt the lexicographic order). The 
box + M(N,H N ) is included in M(2N,H N ) and with high probability b d K (y in) ) 
is connected to a seed in y^ + T(N,h). Nevertheless, piling up the boxes in an 
arbitrary way would give little control on the location of the boxes and therefore 
the connections might not remain in the box M(2N,L' N ). In order to ensure that 
this seed is in T(N, (n+ l)h) one needs to choose a proper subfacet according to the 
steering rule described below. 
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We write y^ = {y\ }i ^ d- Choose j(TT n ) such that the site {—?/[, . . . , —y^Zn h} 
belongs to T j(7Tn) (N, h). By construction + T j(7Tn) (N, h) C T(N, (n + 

where the event C^ n) (N, h) is a translate of the event introduced in Definition 13.31 
and is supported by the domain y( n > + M*(N, H N ). As C 3 K (N,h) is an increasing 
event, we derive a lower bound uniformly over the configurations 7r n in A n 

K"(2N,h n) (4 K) (2iV, hj) > &n* {N>HN) (C&N, h))>\- exp ( - CK*- 1 ) , 

where the last estimate follows from Proposition 13.51 

Proceeding recursively, and applying the steering rule at each step we obtain 

with Ln = L' N i 2 and ip(K, N) = <f'(K, N/2). Using the same strategy as in Propo- 
sition 13. 5[ the position of the seeds can be controlled. This completes the proof. 

□ 

The steering rule was introduced in |BGN| . It will be also extremely useful to 
perform the dynamic renormalization in a two-dimensional slab (see Section 

Remark 4.1. The previous Theorem enables us to generalize inequality (|3.18J) to the 
domain M(N, Ljf). For a given rj > 0, choose K and <p(K, N) according to Theorem 
\4-l\ Then, for any m > 0, there exists M' = M'(rj,m) such that for any N large 
enough, one can find a height ip = <p(r], m, K, N) e [<p(K, N) — M', tp(K, N)] so that 

<S>i* (N , LN) (Y(N, ( p)>m)>l-2r ] . (4.2) 

The derivation of this inequality follows the scheme of the proof of (|3.18|) . 

4.2. Occupied blocks. Theorem 14.11 implies that a connection occurs with high 
probability from a small region around the origin to a height almost at the top of 
a box M(N,Ln). To perform the dynamic renormalization, connections from the 
bottom to the lateral sides of a box will be necessary. This will require different 
arguments. 

We state now the counterpart of Definition 13.31 for the lateral connections. 
Definition 4.1. 

(1) For any integers (£, h), denote by Ck{£, h) the event that 6^(0) is connected by 
open bonds strictly within M(£, h) to a seed lying entirely in S(£, h) (the union 
of the sides). Depending on the orientation of the side, the corresponding seed 
will be of the type fr^-(-) for some appropriate j . 

(2) We also define the collection of events C l K (£, h) for all i ^ (d — l)2 d ^ 1 . Any 
bond configuration in C l K (£, h) is such that b d K (Q) is connected strictly within 
M(£, h) to a seed included in S{£, h) and centered at a site in Si(£, h). 
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(3) A block M(i, h) is said occupied if the bond configuration within this block 
belongs to the intersection of all the sets h) and C 3 K (£, h). 

The main step will be to prove that there is a class of blocks which are occupied 
with high probability. 

Theorem 4.2. Fix i] > 0, then there exists K = K(t]) such that for N large enough 
one can find M(K), h = h(K, N) e [L' N - M(K), L' N ] and I = £(K, N) <: h{K, N) 
so that 

with L' N = L' N (K) ^ 3N. 

The Theorem is a consequence of Propositions 14.11 and 14.21 which are derived in 
the following Subsections. 

4.2.1. Lateral connections. Recall that Y(i, h) denotes the number of sites in T(£, h) 
which are connected strictly within M(£, h) to 6^(0) by open bonds. We also define 
X(£, h) as the number of sites in the sides S(£, h) which are connected strictly within 
h) to b d K {Q) by a path of open bonds. 



Proposition 4.1. Fix 77 > 0, then there exists K = K(i]) such that for N large 
enough one can find (f 1 = <f l {K, N), i 1 = £ X (K, N) and L' N ^ 3N so that 

Vj ^ (d-l)2*-\ ^^(dj^V)) >l-e 1 ( V ), 
where S\{-) converges to as rj tends to 0. 
Proof. 

We start by fixing the parameters according to 77. Remark 14.11 implies that there 
exists K = K(r)) and a sequence L n , such that for any m > one can find M = 
M(r), m, K) and (p n = y? n (?7, m, K) G [L n — M, L n ] such that for any n large enough, 

<*(n,L„) (Y(n, <p n )>m)>l- ^rff- 1 » 1 - r? 2 . (4.3) 
Finally, let m = m(rj, K) be such that (see (|3.19|1 ) 

As m has been chosen large enough, ()4.3|) implies that uniformly over n large 
enough (see Propositions 13.41 and 13. 5j) 

V* < 2 d -\ ®l* {n , Ln ) (C l K (n, ip n )) >1- V 2 . (4.4) 
For any n large enough, there exists £ n such that 

&M*(n,L n ) i Y ^n ~ 1, <fn) ^ ™ ~ l) > 7] , 
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To see that (|4.5|) is well defined, first notice that 

4*(n,L,0 (Y(£ - l,<p n ) ^ m - l) = 1 - & M * {n>Ln) (Y(£ - 1, <p n ) > m) . 

Then it is enough to observe that the function t — > ^L^jyft^n) ^ m ) is 
non-decreasing and 

$ B*(n,L„) ( Y (h <P n ) >rn)=0 and ^. (n>Ln) (Y(n, <p n ) > rnj ^ 1 - rj . 
Thus one defines £ n as the smallest value for which (j4.5j) holds. 

The proof will be decomposed into 3 steps. 

Step 1. 

In this first step, we show that lateral connections occur in a domain smaller than 

B{£ n ,<p n ). 

Lemma 4.1. There exists an integer U = U(i]) such that for any n large enough, 
one can find £' n G [£ n — 1 — U, £ n — 1] and 

where K and (p n were chosen in ()4.3|) and i n was determined in (14. 5(1 . 
Proof. 

We first remark that for £ n defined as above then 

*B.(n,L n )(x(tn-l,<Pn)=0) ^V- (4-6) 

To see this, we apply ()4.3|) and obtain 

V 2 > &m*(n, Ln )(Y(n,<p n ) <m) 

> &m*(n,L n ) (Y(L -l,<p n )<m- 1, X(£ n - 1, <p n ) = 0) 

> &®*(n,L n ) {Y(£n -l,<p n )<m-l) C(n,L„) -l,<p n )=0), 

where the last inequality follows from the fact that both events in the RHS are 
decreasing. Thus inequality ()4.5|) implies ()4.6|) . 

Let us also check that £ n diverges with n. First notice that there exists c p > 
such that 

Vfc < <p n , C(n,L„) {Y(*n, k) = | Y(£ n , k - 1) > 1) > c l ; , (4.7) 

where c e p " is simply the cost of closing all the bonds linking the level h — 1 to the 
level k. 

1 - V < ®l(n,L n ) {Y(tn, Vn) >m)^ 4*(n,Ln) ( Vfc < <Pn, Y(£ n , k) > l) . 

In order to estimate the RHS, we follow the strategy used in Proposition 13.41 We 
proceed recursively by conditioning at each level and using (|4.7|) . 

1 - T) ^ (1 - c^Y" < exp(-Cp n n) , 
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where we used that (f n is 3n. This implies that for n large enough 

L> , 7^7 . log n. (4.8) 
log(l/c p ) 

The derivation of Lemma f4. II follows closely the one of Proposition 13.41 Let u, U 

be two integers which are going to be fixed later. Define recursively the sequence of 
random lengths by Co = £ n — 1 — U and 

C i+1 = inf {/ > d, 1 < X(l, <p n ) < u} A (4 - 1) . 

Notice that for n large enough £ is positive (see ()4.8|) ). 
As in f!3.16|) . we get 

V* < U, $^ {njLn) (A < 4 - l) ^ r^ti)*" 1 , 

where r p (u) < 1 is related to the probabilistic cost of closing at most u open bonds 
(see friSJ)). 

Fix u — u(rj) large enough such that (l — Cp d ^ ( 4K < r//3 (see (|3.19|l ). Then 
choose v = v(r),u) such that r p (u)^ _1 < r//3. Following the derivation of (|3.17|) . we 

get 

$ B*(n,L„) ( Yl hxiLVn) ^ $B*(n,L„) ( X ^n ~ 1, ¥>n) = 0) 

V=*»-l-!7 / 

f 4 

+ $ B*(n,L n ) (£v < in ~ 1) < 3^, 

where we used (14. 6j) . 

Finally as in the derivation of (J3.18|) . one can find U(r), u) large enough such that 
uniformly in n there exists £' n G [t n — 1 — £/, 4 — 1] 

Therefore the probability that a seed lies on top of one of the u attachment sites 
is large. This concludes Lemma 14. 11 

□ 

Lemma 4.2. There exists an integer U = U{rj) such that for any n large enough, 
one can find £' n & [£ n — 1 — U, £ n — 1] and 

Vi < (d - l)2 d -\ & w[niLn) (4(C ^ 1 " (2r/) , (4.9) 

where K and ip n were chosen in ()4.3|) and £ n was determined in (14. 5|) . 

The derivation is similar to the one of Proposition 13.51 

Step 2. At this stage, we know that with high probability 6^(0) is connected to a 
seed in each side of the box 1(4, tp n ) with £' n e [£ n - 1 - U, 4-1]. As 1(4, Vn) 
is contained in M(£ n ,ip n ), we do not have any information on the connections from 
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Figure 4. The events £ ,£ and £ d are depicted. 

6^(0) to the top face T(£' n ,(p n ), thus we need to prove that the lateral connections 
occur as well in a larger box. 

We will use n as the reference scale and consider blocks defined in terms of the 
new parameters (see figure 0} 

t n = [ylogn], Vn = t n + ip n + ipt n , ( 4 -!0) 
N = n + t n + L tn , L' N = ip n + 4ip tn + L tn . 

As (p n ^ 3n (see Theorem 14 .lj) . one has L' N ^ 3N . 

We are going to show that * s connected with high probability strictly within 

M(£' n + 3t n , ipl) to the sides of this block. 

Lemma 4.3. Uniformly over n large enough, 

^l {N , LW) (X(C + 3t n ,^ n ) > 1) (4.11) 

where lim^o 8{v) = 0- 

Proof. Let £ 1 = C x K (t n , ip tn ) be the event that &#(0) is connected to a seed attached in 
T 1 (t n , <y9 tn ) by a path strictly within M(t n , ft n )- From (|4.4jl and the FKG inequality, 
we have 

i(^)(^) >l-V- (4-12) 

For any bond configuration in E 1 , we denote by z^ 1 ' the center of the seed in 
T 1 (i n ,y) 4n ) with the largest coordinate in the lexicographic order. By construction 

Vi ^ d - 1, < zf } < £ n , and zf = <p tn . 
The set S 1 can be partitioned into 

£ l = U ^ (1) )> 

«WeTi(t„, Vtn ) 

where f 1 ^ 1 )) is the set of bond configurations in £ x for which the seed connected 
to 6^(0) with the largest coordinate is centered in z^\ 
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Given z^ in T 1 (t„, (ft n ), we introduce the event S 2 (z^) such that b d K (z<U) is 
connected strictly within z^ + M(£' n , Lp n ) to a seed b^z^) lying in z^ 1 ' + Si(£' n , (f n ) 
with attachment site z^ = {z\ 

4 2) = C + z?> , ^ ^ 4 2) < p tn + V? £ {2, d — 1}, < zf ^t n + £' n . 

Once again, is chosen wrt some arbitrary order and therefore is uniquely deter- 
mined. 

Let be the cr-algebra generated by the bond configurations in M(t n , (p t „). 
Thanks to Lemma f4. 21 we have for any z^\ 

&M*(N,L> n) (£ 2 (^ (1) ) I > <*(n,L n ) ¥>»)) > 1 " , (4-13) 

where we used the FKG property and the fact that for any z^ E T l {t n ,Lp tn ) then 
z {1) + B(n, L n ) remains in B(N, L' N ). 

Similarly, £ 2 (z^) can be partitioned into 

S 2 (z^) = |J S 2 (z^,z^). 

For a given z^ 2 \ let jF z(2> be the cx-algebra generated by the bond configurations in 
the half plane {x, x± ^ zf}. Define the event Ss(z^) that b 1 K (z^) is connected 
to a site z^> which satisfies 

zf - t n < < zf + £ n , Vz e {2, d - 1} , 

(3) 

< ^ ^ (p tn + <fn + tn , 

furthermore the connection must occur strictly within the set {zf, zf + <ft n } x 
{— t n , £' n + 2t n } d ~ 2 x {0, y? n + <f tn + t n }. In the previous computation we used the 
fact that ip n ^ 3n (see Theorem I4.1JI . 

We are going to check that (J4.ll)) is implied by the following inequality 

4£. (W (^V 2) )|^ W ) ( 4 - 14 ) 

By construction 

^M*(N,L' N) (X(£' n + 3t n ,^ n )^l) 

> E E «Ww> n * V 1 *, ^ (2) ) n ^ (2) )) > 

^( 1 )GTi(t n ,^ n ) z ( 2 )e{ z ( 1 )+Si(C^„)} 
where the events in the RHS are disjoint by construction. Applying ()4.14|) we obtain 

<*(n,v n ) + 3t »> rf) ^ X ) ^ E $ b Wjv ) H £ 2 (z (1) )) (1 - rf) . 

Using ()4.13|) and then (|4.12|) . Lemma f4. 31 is satisfied with 

5(rj) = 1 - (l - (277) (d-D^- 1 ) (1 - r?) 2 sC (277) (d-i)^- 1 + 277 . 
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It remains to derive (J4.14j) by piling a block oriented in the direction e\. We 
introduce B«(L, H) = {0, H} x {— L, L} d ~ l and its top face is T«(L, H) = {H} x 
{-L, Ly- 1 . According to l' n > t n for n large enough. Thus ^ (2) + M^(t n , L tn ) 
is included in M(N,L' N ) but does not intersect M(t n ,(pt n ). Inequality (|4.4|) implies 
that with probability larger than 1 — 77, the seed b 1 K (z^) is connected within z^ + 
(t n , (p tn ) to a seed in z {2) + T (1) (t n ,y^J attached at the site z (3) . 

□ 

Step 3. 

Lemma fOl implies that with high probability, b d K {$) is connected to S(£' n +3t n , ipl) 
within M(£' n + 3t n ,cpl). Thus we can proceed as in the first step and derive on a 
larger domain a result similar to (J4.9j) . There exists an integer U = U(rj) such that 
for any n large enough, one can find £\ in [£' n + 3t n — U, £' n + 3t n ] such that 

Vj ^ (d - l)2 d -\ ^(Ar^)(cl(4^n)) ^1-(2^))(^P 3=T . (4.15) 

This completes the proof of Proposition 14.11 with £i(r/) defined by the RHS of the 
previous inequality. 

□ 

4.2.2. Top connections. For any rj, the parameters are fixed according to Proposition 
|4~T1 In particular, K was introduced in (|P|) . iV in (j4~TT)j) and ^ in (|4~T5|) . 

Proposition 4.2. T/jere ezzste £ 2 = £ 2 (K, AT) ancV = ^ 2 (AT, JV) e [L^-M(iT), L' N ] 
such that 

V^2 d ~\ <^ WW) (C]^ 2 V)) ^ 1 - e 2 (77) , 

where £2(0 converges to as rj tends to 0. By construction, £ 2 (K,N) ^ £ Y {K^N) 
and v 2 (K,N) ^ (p l (K,N). 

Proof. Let £ n be defined according to (|4.5|) . Proceeding as in the derivation of (|4.3|) 
and since m has been chosen large enough, we have 

4»(„, in )(Cic(4,^)) ^ 1 - 277 . (4.16) 

We recall that (j4.3|) also implies for n large enough that 

V* < 2 d -\ ^ {tn ^ n) {C K {t n ^ tn )) ^l-rj 2 . (4.17) 

As in (I4.1UJ) . the scaling parameter is n 

£ 2 (K, N)=£ 2 n = £ n + t n , <p 2 (K, N) = V 2 n = 5^ tn + Vn . (4.18) 

According to Theorem 14.11 (p tn ^ 3t n , thus for n large enough 

£ 1 n = £ n -l-2U + 3t n ^ £ 2 n , <pl = t n + <p n + <ft n ^ fl ■ 

Finally, we recall that L' N = ip n +4ip tn +L tn (see fETTUJO so that ^ L' N -ip 2 n ^ M(K). 

In order to connect b d K {Qi) to the height </? 2 , it is enough to link 6^(0) to a seed in 
T(£ n , ip n ) and then to join it to the height y? 2 by piling up 5 blocks of type B(£ n , (ft n )- 
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As in the derivation of Theorem 14.11 the steering rule will be necessary to control 
the position of the seeds. This ensures that the connection from b d K {0) to T(£ 2 , ip 2 ^) 
occurs strictly within B(£ 2 , (/? 2 ). 

Combining inequalities ()4.16|) and ()4.17|) . we derive the Proposition with 

e 2 (v) = 1 - (1 — 277) (1 — 77 2 ) 5 ^ 27] + 5r] 2 . 

□ 

4.2.3. Proof of Theorem^ Setting i(K, N) = ^(K, N) and h(K, N) = <p 2 (K, N), 
the Theorem 14 . 21 follows by combining the FKG inequality and Propositions 14. II and 

5. Slab percolation: second renormalization step 

Using the renormalized blocks, we will prove that percolation occurs in a slab with 
positive probability as soon as p 6 0,. For this, we follow the cluster-growth algo- 
rithm introduced in jBGJNj and prove that on a coarse grained scale large enough, the 
renormalized process dominates a supercritical two-dimensional Bernoulli process. 
This will be enough to ensure percolation in a slab. 

Most of the work was done in the previous Sections to establish uniform estimates 
on the occupied blocks. As a consequence, this final step is very similar to the one 
devised for independent percolation |BGN| . With some respect, it is even simpler 
since the shape of the blocks (ratio height /width) is under control. 

Let us fix p G Q q and summarize the result of Theorem 14.21 For any 77 > 0, there 
exists K, (L, H) and (£, h) such that 

®]&*(L,H) (®*(^> h) is occupied) ^ 1 - rj , (5.1) 
where K characterizes the size of the seeds and the parameters are chosen such that 

H^3L, H- h ^ tf/100. 
At this stage the dependency on N is no longer relevant. 

We sketch below the main steps of the proof and detail only the new features. 

5.1. Reduction to two dimensions. The strategy will be to create new connec- 
tions by stacking rotated translates of the block B(£, h). In the derivation of Theorem 
14.11 we already encountered two basic features: top stacking and steering. We 
will also present a third one: branching. This will enables us to define a new 
dependent percolation process restricted to the slab Sl — {— 2L, . . . , 2L} d ~ 2 x Z 2 . 

Steering is necessary to control the deviation of the renormalized paths. By choos- 
ing carefully on which particular subfacet of a block another block is attached, we 
may localize a trajectory. In particular, the first (d — 2) coordinates happen to be 
irrelevant because any stacking in a direction el (for i ^ d — 2) can be centered along 
the i-axis (we refer to [BGN for a complete explanation). We stress the fact that by 
steering the sequence of boxes ~B(£, h) will remain inside = {— 21, ■ ■ ■ , 2£} d ~ 2 x Z 2 , 
nevertheless to evaluate the probability of an occupied block one needs to average 
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over a bigger block M(L,H) (see (j5.1j) ). This explain why we have to consider the 
renormalized process in the thicker slab Si+e = {—L — £,..., L + £} d ~ 2 x Z 2 . 

Let us now introduce some notation to emphasize the role of the plane (e^-i, e^) 
which we shall refer later on as the (x, y)-plane. In order to describe the geometrical 
construction in the (x, y)-plane, it is convenient to rename the subfacets of the 
blocks. The direction will be referred as North, — as South, Cd-\ as East and 
— e*rf_i as West. In particular M(£, h) will be dubbed north block and denoted by 
B North (£, h). Its north face T North (£, h) = T(£, h) is split into 

T N orth,East = {-1,1}*'* * {M} >< W , 
T N orth,We S t = {-£,£} d ~ 2 X {-£,0} X {H} . 

We also distinguish the western and eastern sides 

^North.East = {~£,£} d - 2 X {£} X {0,#}, 
^No rt h,West = {-£,£} d ~ 2 X {-£} X {0,H}. 

Rotating B North leads to define blocks oriented in the other 3 directions B East , Bg outh 
and Bwest- Their faces will be named by transposing the previous notation. 

We can now describe the main ingredients to concatenate paths from different 
occupied blocks. Starting with an occupied north block, top stacking enables to 
pile up n occupied blocks in the north direction with a probability of success at least 
(1 — rj) n . As this procedure is doomed to fail sooner or later, a branching procedure 
will be necessary to allow percolation in the other directions. 

An occupied block contains attachment sites on each of its faces, so that several 
blocks can be stacked on it. A branching occurs when two blocks are stacked 
simultaneously on the top and on one side of an occupied block. These blocks should 
be positioned in a careful way to remain essentially independent. For example, 
starting with the occupied block B North (£, h), a block y + B North (£, h) can be stacked 
on a seed centered in y G T NoTthjEast (£, h) and a west block 2 + B West (£, h) on a seed in 
5 , North,West ( see figure EJ). Alternatively, exchanging East and West would lead to a 
branching in the North/East directions. By construction, the event that these new 
blocks are occupied, is supported by the disjoint set of bonds y + B^ orth (£, h) and 
z + B^g^f, h). Nevertheless to evaluate the probability of a block to be occupied 
requires averaging over a larger domain M,(L,H) (see (j5.1|) ). This raises some 
measurability issues which will be detailed in the example below. 

We are going to evaluate the probability that the following sequence of occupied 
blocks occurs (see figure EJ). The derivation will be similar to the one of Lemma f4. 31 
In the procedure below, the steering in the first (d — 2) coordinates is implicit. 

First 6^(0) is connected strictly within BN or th(^, h) to a seed centered at the site 
which belongs to T^ 0Tt h^ est (£,h). The site is uniquely determined and is 
chosen as the first site (wrt to the lexicographic order) for which the previous event 
holds. This event is denoted by £ 1 (y^). Thanks to the key estimate (|5.1|) . we have 
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Figure 5. A branching sequence in the North/West directions. The domains 
+ B North (L, H) and z« + M West (L, H) are depicted by the dashed blocks. 

for any domain A containing M(L, H) 

$a (£\y (1) )) = $A (Ck& h)) > &n* {L , H ) (B*(4 fc) is occupied) ^ 1 - 77 , (5.2) 

where the sum over t/ 1 ) is restricted to T NorthiWe st(^ 

Then a branching occurs within y^ 1 ' + B North (£, h), i.e. bj < (y^) is connected to a 
seed bx(y^) centered in y^ 2 > £ y^ 1 ' + T North)East (£, h) and b d K (y^) is also connected 
to a seed b d ^ 1 (z^) centered in G y^ + S^ OTt h,west(^, h) . This event is denoted 
by E 2 (y( l \ y^ 2 \ z^). As the attachment sites y^ 2 \ z^ 1 ' are uniquely determined, the 
previous events are disjoint for different attachment sites. 

Finally, the connections evolve in the West and North directions: b d ^ 1 (z^) is 
connected to a seed b d ^ 1 (z^) in z^ + Tw est (l, h) and bj ( (y^ 2 ') is connected to a seed 
in y (2) + T NorthjEast (£, h) or possibly to a seed in y (2) + SWorth.East^, h) (on figure El 
a side stacking from ,2® oriented in the East direction is also depicted). We also 
define 

£3 (y (2)) = | y ( 2 ) + B North (£ ) ft) i s occupied} , 
+ Bwest(^) h) is occupied} . 
We are going to check that a branching occurs with high probability 

J2 *i(^(y (1) ) n^^w,^) nfV) n ^W)) ^ (1 - „)*. (5.3) 

It is enough to iterate in the right order the key estimate (|5.1J1 . Let us fix a triplet 
y( 2 ) ? ^C 1 ) } and drop temporarily the dependency on the sites in the events 
£'. The event £ l fl £ 2 fl £ 3 is supported by bond configurations in the hyperplane 
{%d-i ^ Vd-i — and the support of £ 4 lies in {xd-i < y^\ — £}■ Thus conditioning 
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outside + B West (L, H) and using the fact that £ 4 is non decreasing, we get 

> $1 (£ x n £ 2 n £ 3 ) $ f ,( 1)+Bwest(ii H) n > $a n £ 2 n £ 3 ) (i- v ). 

At this stage there are no more ambiguities and we estimate the remaining events 
one after the other (starting from the top) thanks to (|5.1|) . This completes (15. 

5.2. Cluster-growth algorithm. We describe now the second level of renormal- 
ization. 

As explained in the previous Subsection, it is enough to consider a two-dimensional 
projection of the system. The projection of the block B.(£, h) onto the plane (x,y) 
will be called a brick and denoted by B.(£, h). The orientation convention applies 
as well for the bricks. To an occupied block is associated a successful brick with 
four connection sites one in each set T North)East , T Nort h,West, ^North^ast, ^North.West 
which are the projections of the centers of the seeds. A sequence of bricks build 
by concatenating occupied blocks according to the stacking procedures previously 
described will be called a successful brick sequence. 

Thus there is an immediate correspondence between blocks and bricks and it is 
enough to build the second level of renormalization in terms of bricks in Z 2 . The 
second renormalization level lives on the coarse grained scale iV = 10L + 10H. 
The lattice Z 2 is partitioned into translates of the square So = {—N + 1, . . . , N} 2 
which are indexed wrt an arbitrary order {§j}j ^o- The algorithm is performed by 
inspecting the squares one after the other and checking iteratively some properties 
which will be detailed later on. If these properties are satisfied then the square 
contains a crossing cluster made of successful bricks. A random variable Zj(a>) 
depending on the bond configuration u is associated to the square Sj. The square S; 
is declared to be good if these properties hold and we set Zj = 1, otherwise Zj = 0. 
As it will be clear later on, any bond configuration must contain an open cluster 
intersecting all the good squares. 

We explain now the construction of the renormalized process {Zj}j indexed by 
the squares in the domain {— M, . . . , M} 2 for some M = 2mN. Let us suppose that 
k squares So, S^, . . . , Sj fc _ 1 have been examined. Choose the next square E>i k as the 
earliest square (in the fixed ordering) which has a face in common with one of the 
previous good squares. If no such square exists then the algorithm stops, otherwise 
we are going to prove that for a suitable tuning of the coarse grained scales 

*5 w (^ fc = 1 l^*-i'---^o)^a ) (5.4) 

where a = a{rf) can be chosen arbitrarily close to 1 and in particular much larger 
than the critical value of the site percolation in Z 2 . We stress the fact that a depends 
only on i], so that all the parameters are determined thanks to (|5.1|) . Inequality 
(J5.4J) implies that the cluster formed by the good squares dominates stochastically 
the open cluster growing from the origin of a supercritical site percolation process. 
We refer to Grimmett, Marstrand jGMj for a detailed account of the stochastic 
domination. 
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Before completing the construction of the renormalized process and deriving (J5.4j) . 
we pause to look at the proof of Theorem 12.11 

Proof of Theorem \2.1\ 

Fix M = 2mN and choose a site x in iSl+^m- Let S be the square such that the 
tube T x = {—L - £, L + £} d ~ 2 x § contains x. If x = 0, then T = {—L -£,L + 
£} d ~ 2 x So- Define Z(0,x) as the set of bond configurations such that there exists 
a renormalized connected path of good squares from S to S. By construction, one 
has 

This means that the tubes are connected by an open path if a connection occurs on 
the coarse grained level. Thanks to the stochastic domination (see (J5.4j) ). one can 
choose a = a(rj) such that the RHS is bounded from below uniformly wrt to x and 
M. Therefore p ^ p c (L +£) ^ p c and Theorem 12. II holds. 

□ 

We turn now to the derivation of ()5.4|) . As explained previously, the cluster- 
growth algorithm examines in turn each square: a square § will be inspected only if 
there exists a square §' which shares a common facet with S and has been previously 
declared to be good. In this case the state of S will be determined independently 
of the bonds lying outside §U §' (more precisely, only a small portion of S' will 
be relevant). The procedure we are going to apply is translational and rotational 
invariant, thus it is sufficient to describe it in a special case: conditionally to the 
fact that §o is a good square, we would like to determine the state of its northern 
neighbor Si = S + (0, N). 

We define the target regions of So as 

T° OIth = {-N/2,N/2}x{N-2H} 

and the other target regions ^ ast , T^ outh , T^ est are deduced by rotation. In the same 
way, the target regions of Si are obtained by translation Tj- = + (0, N). We will 
see that a necessary condition for a square to be good is to contain a successful brick 
M m (£,h) intersecting the corresponding target region T,, where • ranges over the 4 
directions. 

Conditionally to the fact that So is a good square, there exists by definition a 
successful north brick B intersecting the target region T^ OTth . This brick will be used 
as a starting point to launch connections to the target regions of Si. The main 
ingredients are the centering and the bifurcation rules (see Figure 

Centering rules: 

First, bricks are piled up on top of IB in the north direction using the steering rule 
in order to center the brick sequence along the axis x — 0. At some point a brick, 
denoted by B will intersect the level {y = N + N/2}, this triggers the bifurcation. 
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triggering line 



T 1 



East 



y 



X 



Figure 6. A successful sequence of bricks is depicted starting from the northern 
part of § and invading Si. The dashed lines are the target regions. 

Bifurcation: 

On top of B lies a connection site = (xi,yi) (recall that z^ is simply the 
image of the center of the seed in the top face of the occupied block associated to 
B ). By construction Xi belongs to {—5H — 5L, 5H + 5L} and we also suppose that 
X\ < 0, the other case can be treated by symmetry. 

First a branching occurs in the North/West directions, to do so three bricks are 
piled up in the north direction. 



From B 2 a brick sequence branches in the West direction towards the target region 
7^ est . We will come back to this later. 

Then an east brick B 4 is attached to the east side of B 3 and another branching 
occurs in the East/North directions: three bricks are piled up in the East direction 



Bi = zV +M North (£, h), 

B 2 = Z& + B North (£, h) With Z& G Z W + T No rth,West(^ h) , 
B 3 = + B North (£, h) With Z^ e + T N orth,East(^, h) . 



on top of B 4 



B 4 = Z^ + B East (£, h) With G S N o rt h,East + , 

B 5 = z^ + B East (£, h) with z^ G T East , North + , 

B 6 = z^ + M East (£, h) with z^ G T East ,No rt h + , 

B 7 = + B East (£, h) with z^ e T EastiSouth + z^ . 



From B 6 a brick sequence branches in the North direction towards the target region 
7^ orth and another sequence goes on in the East direction towards T^ast- 
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Final connections: 

After the bifurcation, the steering rules are applied once again to center along the 
axis of the target region 7^ est (resp. North, East) the brick sequence piled up in the 
West direction (resp. North, East directions). It remains to check that using this 
construction the sequences reach the target regions and that they do not overlap. 

The starting point of the western sequence z^ w ' = (xw, Vw) belongs to the western 
side of B 3 . By construction 

x w G [-21, t\ + x 0} y w e[N + N/2 + H,N + N/2 + 3H] 

By piling up bricks in the western direction, the brick sequence is localized in the 
strip {— N, xw} x {yw ~ L,L} and therefore it encounters the target region %^ est 
by using less than 20 bricks. 

We proceed in the same way for the other directions and see that the sequences 
evolve in non overlapping regions. 

The inspection of a square requires less than 100 bricks, thus the probability 
for Si to be good conditionally to the fact that §o is good, is larger than a{rf) = 
(1 — r/) 100 . Furthermore, the procedure described above does not depend on the bond 
configurations outside SoflSi. Thus the stochastic domination inequality ()5.4|) holds 
at any step of the cluster algorithm. 

References 

[Al] M. Aizenman, Rigorous studies of critical behavior. II. Statistical physics and dynamical sys- 
tems, Progr. Phys., 10, Birkhauser Boston, Boston, MA, 453-481 (1985). 

[A2] M. Aizenman, Absence of an intermediate phase for a general class of one-component ferro- 
magnetic models, Phys. Rev. Lett. 54, no 8, 839-842 (1985). 

[ACCFR] M. Aizenman, J.T. Chayes, L. Chayes, J. Frohlich, L. Russo, On a sharp transition 
from area law to perimeter law in a system of random surfaces, Comm. Math. Phys. 92, no. 
1, 19-69 (1983). 

[BGN] D. Barsky, G. Grimmett, C. Newman, Percolation in half-spaces: equality of critical den- 
sities and continuity of the percolation probability, Probab. Theory Related Fields 90, no. 1, 
111-148 (1991). 

[B] T. Bodineau, The Wulff construction in three and more dimensions, Comm. Math. Phys. 207, 
197-229 (1999). 

[BIV1] T. Bodineau, D. Ioffe, Y. Velenik, Rigorous probabilistic analysis of equilibrium crystal 

shapes, J. Math. Phys. 41, No.3, 1033-1098 (2000). 
[BKM] C. Borgs, R. Kotecky, I. Medved, Finite-size effects for the Potts model with weak boundary 

conditions, J. Statist. Phys. 109, no. 1-2, 67-131 (2002). 

[C] R. Cerf, Large deviations for three dimensional supercritical percolation, Asterisque No. 267 
(2000). 

[CP] R. Cerf, A. Pisztora, On the Wulff crystal in the Ising model, Ann. Probab. 28, no. 3, 
947-1017 (2000). 

[DS] R. L. Dobrushin, S. Shlosman, Completely analytical Gibbs fields, Progr. Phys., 10, 
Birkhauser Boston, Boston, MA, 371-403 (1985). 

[ES] R. Edwards, A. Sokal, Generalization of the Fortuin-Kasteleyn-Swendsen- Wang representa- 
tion and Monte Carlo algorithm, Phys. Rev. D (3) 38, no. 6, 2009-2012 (1988). 

[FK] C. Fortuin, P. Kasteleyn, On the random- cluster model. I. Introduction and relation to other 
models, Physica 57, 536-564 (1972). 



SLAB PERCOLATION FOR THE ISING MODEL 



33 



[FP] J. Frohlich, C. Pfister, Semi-infinite Ising model. II. The wetting and layering transitions, 

Comm. Math. Phys. 112, no. 1, 51-74 (1987). 
[Gl] G. Grimmett, Percolation, Second edition, 321, Springer- Verlag, Berlin, 1999. 
[G2] G. Grimmett, The stochastic random- cluster process and the uniqueness of random- cluster 

measures, Ann. Probab. 23, no. 4, 1461-1510 (1995). 
[GH] G. Grimmett, P. Hicmcr, Directed percolation and random walk, Progr. Probab., 51, 

Birkhauser Boston, Boston, MA, 273-297 (2002). 
[GM] G. Grimmett, J. Marstrand, The supercritical phase of percolation is well behaved, Proc. 

Roy. Soc. London Scr. A 430, no. 1879, 439-457 (1990). 
[KLMR] R. Kotecky, L. Laanait, A. Messager, J. Ruiz, The q-state Potts model in the standard 

Pirogov-Sinai theory: surface tensions and Wilson loops, J. Statist. Phys. 58, no. 1-2, 199-248 

(1990). 

[LP] J. Lebowitz, C. Pfister, Surface tension and phase coexistence, Phys. Rev. Lett. 46, no. 15, 
1031-1033 (1981). 

[MS] R. Meester, J. Steif, On the continuity of the critical value for long range percolation in the 

exponential case, Comm. Math. Phys. 180, no. 2, 483-504 (1996). 
[MMP] A. Messager, S. Miracle-Sole, C. Pfister, On classical ferromagnets with a complex external 

field, J. Statist. Phys. 34, no. 1-2, 279-286 (1984). 
[MMR] A. Messager, S. Miracle-Sole, J. Ruiz, Surface tension, step free energy and facets in the 

equilibrium crystal, J. Stat. Phys. 79 (1995). 
[Pi] A. Pisztora, Surface order large deviations of Ising, Potts and percolation models, Prob. Th. 

Rel. Fields 104, 427-466 (1996). 
[SS] R. Schonmann, S. Shlosman, Complete analyticity for 2D Ising completed, Comm. Math. 

Phys. 170, no. 2, 453-482 (1995). 

Departement de mathematiques, Universite Paris 7, case 7012, 2 place Jussieu, 
Paris 75251, France 

E-mail address: bodineau@gauss.math.jussieu.fr 



